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Abstract 

The behaviour of semi-dilute wormhke micelle solutions under vertical vibrations is 
investigated using classical measurements of the Faraday instability threshold along 
with two new experiments based on particle imaging velocimetry and birefringence. 
We provide evidence for the presence of oscillations of the critical acceleration and 
wave number with the vibration frequency, which are linked to oscillations of the 
velocity and birefringence fields along the vertical direction. These observations 
are interpreted in terms of standing elastic waves across the fluid layer. Such an 
interpretation is confirmed numerically by using the model proposed by S. Kumar 
[Phys. Rev. E, 65, 026305 (2002)] for a viscoelastic fluid to calculate the velocity and 
deformation fields. Finally, further birefringence experiments above the instability 
threshold show that the Faraday instability does not induce disentanglement and 
breakage of the micellar network. 

Key words: Wormlike micelles, Faraday instability, particle imaging velocimetry, 
birefringence 



Introduction 



In the classical "Faraday experiment" a fluid layer is submitted to vertical 
vibrations at a frequency / and acceleration a ([H). By increasing the accelera- 
tion from zero, the system undergoes a bifurcation from a simple, well-defined 
state (fiat interface with zero velocity in the bulk) to a more complex, desta- 
bilized state: at a given acceleration Oc called the critical acceleration, the 
fiat interface goes unstable and a surface wave pattern characterized by the 
critical wave number kc develops ([H; 0; 0; 0). When the driving acceleration 
is further increased above Oc, secondary instabilities take place (jsl), leading 
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at even higher accelerations to spatio-temporal chaos (u). Together with the 
Rayleigh-Benard convective instability, the Faraday instability has often been 
described as a paradigm for the study of dynamical systems at the laboratory 
scale (lZi). 



Until about a decade ago, the Faraday instability had been intensively studied 
in Newtonian fluids (0; HI; 0; H; 0) . Only recently have vertically vibrated com- 
plex fluids been the subject of experimental (llOl : I 111 : Il2l : Il3l : Il4j ) and theoretical 
research Hot . \lw . It was shown that the coupling between the microstructure 
of complex fluids and the Faraday instability either only slightly modifies the 



instability (e.g. by affecting the critical acceleration (llOl : I 111 ) or by delaying 



the transition to disordered states (Il4l)) or presents non-intuitive behaviours 
such as stabilized "persistent holes" in shear-thickening fluids (1l2l). With the 
notable exception of concentrated colloidal and granular suspensions (|l^: [itI). 
most of these works were performed on fluids whose viscoelasticity is only a 
small perturbation (fiol : [isf : 14) to Newtonian behaviour. The case of a com- 
plex fluid with significant viscoelasticity was addressed in Ref. (jlli ) and it 
was shown that viscoelasticity could lead to a harmonic response of the fluid 
surface {i.e. at /) instead of the classical subharmonic response {i.e. at //2). 



In a previous paper (llSl ) we reported onset measurements of the Faraday in- 
stability in a semi-dilute surfactant solution (CPyCl-NaSal at 4% wt.) known 
to form an entangled network of "wormlike" micelles. In this very viscoelas- 
tic fluid, the critical acceleration and wave number were shown to present 
oscillations as a function of the driving frequency /. We interpreted these os- 
cillations in terms of standing elastic waves between the disturbed surface and 
the container bottom. 



The present article is devoted to a more thorough study of this striking ef- 
fect of viscoelasticity on the Faraday instability. After a brief description of 
our wormlike micellar system and of our experimental setup, we flrst present 
new sets of Faraday experiments in various wormlike micelle solutions and 
explore the dependence of the fluid response on surfactant concentration and 
on temperature. Then, the interpretation in terms of standing elastic waves 
is conflrmed using two novel experiments: particle imaging velocimetry (PIV) 
and birefringence under vertical vibrations, which show a non monotonic be- 
haviour of the velocity and deformation flelds along the vertical direction. 
We also use the flnite-depth model of ref. (llSl ) to predict the velocity fleld. 
This numerical calculation accounts qualitatively for all our experimental ob- 
servations. Finally, reasons for quantitative discrepancies between the model 
and the experiments are discussed. In particular, birefringence measurements 
above onset show that the speciflc nonlinear features of the rheology of worm- 
like micelles, namely disentanglement and breakage of the micellar network 
and shear-induced alignment, do not come into play in our experiments. 
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1 Viscoelastic fluid: wormlike micelle solution 



Wormlike micelles are long (typically hundreds of nanometers to tens of mi- 
crons), cylindrical (a few nanometers in diameter) and flexible aggregates of 
surfactant molecules in aqueous solution (jl9l ). Depending on the concentration 
and on the surfactant, wormlike micelle solutions present various character- 
istic behaviours: shear-thickenning at low concentration (l20l) . shear-induced 
alignment in the semi-dilute regime (1211 ). and equilibrium isotropic to nematic 
phase transition at even larg er concentrations (1221 ). In the following we focus 
on semi-dilute solutions (1231 ) . As shown by Gates (124), the unique viscoelas- 
tic behaviour, namely almost a perfect Maxwell fluid (1251 ). of these solutions 
is due to the interplay of two processes: reptation motion as in conventional 
polymers (1261 ) and breaking and recombination of the micelles under thermal 
agitation. Therefore, the complex viscosity of these solutions can be written: 



1 + iuJT2 

Vo-r—. , 

1 -|- tUOTi 



where rjo is the fluid zero-shear viscosity, ri the relaxation time, uj the pulsa- 
tion, and T2 a second characteristic time. 

For a pure Maxwell fluid, T2 = 0. However, due to the presence of the solvent, 
we shall rather use T2 = ti 775/770, where rjs -C rjo is the solvent viscosity (in our 
experiments, rjs = rjbrine = 10~^ Pa.s and rjo ~ 10 Pa.s). Thus, in the case of 
a Maxwell fluid, there are really only two free parameters t/q and ti in eq. (1) 
since T2 simply derives from the solvent viscosity 7/5. Another widely used 
model for the linear rheology of wormlike micelles is the Oldroyd fluid whose 
complex viscosity is also given by eq. (1) but where T2 is left as a free parameter 



(see e.g. refs. ( 1271 : |28|)). We shall see below that this additional free parameter 
allows one to account for the behaviour of our solutions over a much larger 
frequency range than the Maxwell fluid. Note that such a use of the Oldroyd 
model assumes the validity of the Cox-Merz rule which may be violated in 



some wormlike micelle solutions for ujti > 1 (1271 : l29l ). We still use fits by 
the Oldroyd model since equation (1) yields good empirical representation of 
our data. Also note that more elaborate models that account not only for 
the solvent but also for high frequency relaxation modes based on Rouse or 
Zimm approaches may apply to wormlike micelle solutions (^(i). Finally we 
introduce the shear modulus as the high-frequency limit of the imaginary part 
of -ujrj*{Lj) I.e. Go = lim^^oo 3(-a;?7*(a;)) = 7/o(ri - T2)/t^ ~ Vo/n. 

In the following, we study various wormlike micelle solutions made of cetylpyri- 
dinium chloride (CPyCl, from Aldrich) and sodium salicylate (NaSal, from 
Acros Organics) dissolved in brine (0.5 M NaCl) with a fixed concentration 



ratio [NaSal]/[CPyCl]=0.5 as described in refs. (l2l|; M). The CPyCl-NaSal 
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c (% wt.) 


2 


3 


4 


6 


8 


m (Pa.s) 


2.97 


12.8 


27.6 


74.1 


129 


n (s) 


0.421 


0.713 


0.882 


1.08 


1.11 



Table 1 

Best fit parameters for a Maxwell fluid (eq. (1) where T2 = Tirjs/rjQ is fixed) for 
CPyCl-NaSal solutions of various concentrations c. 



c (% wt.) 


2 


3 


4 


6 


8 


m (Pa.s) 


2.95 


12.6 


27.7 


77.5 


129 


n (s) 


0.411 


0.714 


0.891 


1.13 


1.10 


r2 (ms) 


2.00 


0.754 


0.611 


0.344 


0.849 



Table 2 

Best fit parameters for an Oldroyd fluid (eq. (1) where both ri and T2 are free 
parameters) for CPyCl-NaSal solutions of various concentrations c. 

concentration of our solutions varies from 2% to 8% wt. This concentration 
range ensures that (i) the micelles are in the semi-dilute regime, i.e. remain 
entangled and form a viscoelastic network (1211 ) and (ii) the solutions are not 
too viscous (r/o — GqTi < 100 Pa.s) so that the onset of Faraday waves can be 
observed. 



Linear rheological measurements were performed in the cone-and-plate geom- 
etry using a shear rate controlled rheometer (ARES, TA Instruments) with 
small shear rate oscillations of amplitude 0.1 s~^ (the linear regime extends 
up to at least 1 s~^ for all concentrations under study). Such measurements 
yield reliable values for the viscoelastic moduli G'{uj) = —u!Q[ri*{u))] and 
G"{u!) = to'3?[?7*(to')] as long as the oscillation frequency / = uj/2'71 remains 
smaller than about 40 Hz. Figure 1 shows the dimensionless data G' /Gq and 
G" / Go as a function of ujti for the various solutions used in the present work. 
The fact that the G"/ Gq vs uti data do not collapse on a single curve at high 
frequency provides a clear indication that the Maxwell model fails for ujti > 3. 
Indeed both G' and G" are well accounted for by eq. (1) over the whole fre- 
quency range provided T2 is left as a free parameter (Oldroyd fluid). Tables 1 
and 2 gather the best fit parameters using eq. (1) as a function of surfactant 
concentration. As shown by the dashed lines in Fig. 1 (see also Fig. 5(a) for 
raw data measured on a 3% wt. CPyCl-NaSal solution and the corresponding 
Maxwell and Oldroyd fits), the Maxwell model (where T2 = Tirjs/rio is fixed) 
leads to a good description of the G"/ Gq data only for uti smaller than about 
3. This means that higher frequency modes, which are captured by the Ol- 
droyd model, come into play for uti > 3, a feature that is classically observed 
in wormlike micelles (1301). Finally a power-law fit of our t]o vs c data yields 
an exponent of 2.7±0.3 in rough agreement with the value of 3.3 reported 
in refs. (1231 : |22| ) and consistently with the "fast-breaking" limit for micellar 
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Fig. 1. Linear rheology of 2, 3, 4, 6, and 8% wt. CPyCl-NaSal solutions, (a) Storage 
modulus G'/Gq and (b) loss modulus G" /Gq versus wri, where uj is the shear 
oscillation pulsation, ri the G'{uj) the 

shear modulus. The dashed line and the continuous line respectively represent a 
Maxwell fluid (?7o = 49 Pa.s, ti = 0.84 s) and an Oldroyd fluid (t^q = 50 Pa.s, 
n = 0.85 s, T2 = 0.91 ms). 



dynamics (1241 ) 



2 Critical acceleration and wave number measurements 



2.1 Experimental setup 



Our experimental setup is similar to the one used in refs. (js]; 18). It consists 
of a cylindrical Plexiglas container of diameter o? = 60 mm filled to a height 
h = mm under the brimful boundary condition This cell is sealed by a 
Plexiglas cover that prevents evaporation and surface contamination, and ver- 
tically vibrated by an electromagnetic shaker (Ling Dynamic Systems V406). 
Except when stated otherwise, the fluid temperature is controlled to 21±0.5°C 
by water circulation beneath the container. The acceleration is measured by 
a small piezoelectric accelerometer (Endevco 2224C) attached to the cell. The 
wave number is inferred from images of the surface taken by a CCD camera 
(Cohu). The illumination technique allows us to detect surface deformation 
down to an amplitude of about 30 /xm (l3ll ). We checked that our container 
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is large enough for the system to be considered as laterally unbounded by 
performing experiments in various Newtonian fluids (silicon oils and water- 
glycerol mixtures) and by comparing the experimental measurements of Qc and 
kc to the numerical calculation of ref. (@) for an unbounded viscous fluid. As 
already mentioned in ref. (fisl ). the agreement between experiment and theory 
is almost perfect over the whole range of frequencies / = 20-120 Hz, so that 
finite-size effects due to lateral boundaries are negligible. 



2.2 Measurements of and k^. 



At a given driving frequency /, the instability threshold is determined as 
follows: (i) the acceleration is quickly increased until the surface is fully desta- 
bilized then quickly decreased until the instability completely disappears; (ii) 
from this last value, the acceleration is slowly increased again (by approx. 
1%/min) until the interface goes unstable, which yields an "upper limit" acmax 
for the critical acceleration; (iii) once the instability appears, we wait for the 
whole surface to be fully destabilized (which takes about 1 min); (iv) a picture 
of the surface is taken from which the critical wave number kc is estimated; 
(v) the acceleration is finally slowly decreased (by approx. 1%/min) until the 
instability completely disappears, which yields a "lower limit" acmin for the 
critical acceleration. In all cases, the upper and lower limits of Oc differ by less 
than 2%: no significant hysteresis is observed and we define as the average 
o-c = {cLcmin + '^cmax) l'^- Fo^ ^11 our solutious and over the whole range of in- 
vestigated frequencies / = 20-120 Hz, the surface response was found to be 
suhharmonic i.e. the standing wave pattern oscillates at //2. 



As reported in ref. (18:), the critical acceleration and wave number in worm- 
like micelles strikingly differ from their equivalent in Newtonian fluids. Indeed 
both curves and kc present pronounced oscillations as a function of the 
driving frequency /. Figure 2 shows Oc and kc vs / for a 3% wt. CPyCl- 
NaSal solution (see also fig. 1 of ref. (jlSl ) for similar measurements on a 
4% wt. solution). A consequence of this peculiar behaviour is that kc does 
not follow the usual dispersion relation for gravito-capillary surface waves: 
oj"^ = {gkc + crfc^/p) tanh(fcc/i), where u = nf for subharmonic response, g 
is the acceleration due to gravity, a the surface tension and p the fluid den- 
sity (0). Instead we find a non-monotone dispersion relation. Since the same 
experiments performed on Newtonian fluids with similar viscosity and surface 
tension do not display such behaviour, the oscillations observed in fig. 2 must 
be linked to the viscoelastic properties of our micellar solutions. 
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Fig. 2. (a) Critical acceleration Oc and (b) critical wave number kc plotted against 
the vibration frequency / for a 3% wt. CPyCl-NaSal solution at T = 21°C. Error 
bars on Oc are estimated from the experimental "upper" and "lower" limits for ac by 
±{acmax~0'cmin)/'^ (scc text). Error bars on kc correspond to the standard deviation 
of 3 to 6 estimations performed on various areas of the picture of the surface. 

2.3 Simple physical interpretation of the oscillations 




z 

.• . J" 



In ref. ( llSl ) we proposed to interpret the oscillations in ac and kc as the signa- 
ture of "standing elastic waves" across the cell. Indeed in a strongly viscoelastic 
fluid, the disturbed surface is likely to generate shear waves that propagate 
downward across the fluid layer, are reflected at the bottom of the cell, back- 
propagate toward the surface where they may interact with the surface wave 
provided the attenuation over a distance 2h is not too large. The velocity of 
shear waves is given by c{u) = \JG'{uj)/ p ~ ^Go/p for a Maxwell or an 01- 
droyd fluid in the ujti ^ 1 limit, which is always verified in our experiments 
(from table 2, ujti = 27r/ri is in the range 25-400). Using this simplification 
one can calculate the phase of the reflected shear wave as it reaches the surface. 
Indeed, the total travel time from the surface to the bottom of the cell and 
back to the surface is 2h/c, which corresponds to a phase shift (2/i/c)(a;/2) 
since the response of the surface is subharmonic. Moreover the no slip bound- 
ary condition at the bottom of the cell leads to an additional phase shift of vr so 
that the phase of the reflected shear wave at the surface is (p = uh^p/Go + vr. 
li if = 2-/™, with n an integer, the reflected shear wave and the surface wave 
are in phase leading to constructive interference, amplification of the surface 
wave, and a local minimum of ac- On the other hand, when (p = (2n + l)7r, the 
interference is destructive and surface waves are hindered, leading to a local 
maximum of Oc, hence the oscillations of Fig. 2(a). 
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Fig. 3. (a) Critical acceleration Oc versus f /5fth with 6fth = \JGqI pjh for a 3% wt. 
(•) and a 4% wt. (o) CPyCl-NaSal solution at T = 21°C. (b) Critical acceleration 
Oc versus / for a 4% wt. CPyCl-NaSal solution at 23°C (A) and 19°C (v)- 

This simple interpretation, which only focuses on a single reflection and ne- 
glects attenuation of the shear wave, leads to the following expression for the 
distance 5/ between two maxima of a^- 



th 




(2) 



The experiments presented in ref. ( 1l8l ) were restricted to only one sample at 



a single surfactant concentration of 4% wt. The main objective of the present 
work is to further check the prediction of equation (2) over the whole semi- 
dilute regime. More experiments were thus performed by varying the surfac- 
tant concentration, i.e. for different elastic moduli Gq. Figure 3(a) shows the 
critical acceleration versus the reduced frequency f /Sfth for 3% and 4% wt. 
CPyCl-NaSal solutions, where the value of Gq in eq. (2) is given by linear rhe- 
ological measurements. In this representation, critical accelerations oscillate 
roughly with the same period as a function of reduced frequency. However the 
maxima of ac do not exactly correspond to integer multiples of f/Sfth, which 
points to a first limitation of the simple physical picture proposed above. We 
shall come back to this discrepancy in sect. 3.2. 

Equation (2) also predicts that 6f does not depend on the fluid relaxation 
time Ti. For CPyCl-NaSal wormlike micelle solutions, it is well known that 
Ti greatly depends on the temperature while Go is nearly temperature in- 



dependent (1301 ). Indeed on our 4% wt. mi cellar solution, linear rheo logical 



experiments gave Gq = 33 Pa and Ti = 0.74 s at T = 19°C (smaller tempera- 
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Fig. 4. Experimental 5f (•) versus predicted 6fth = \/GqI pjh. The continuous line 
is the first diagonal. 

tures led to crystallization of the sample), and Gq = 32 Pa and ri = 0.37 s at 
T = 23°C: a 4°C change in temperature only induces a 3% change in Gq while 
Ti is doubled. Figure 3(b) shows Oc measurements performed at T = 19°C and 
T = 23°C (slightly below room temperature to prevent condensation on the 
Plexiglas cover). As expected from eq. (2) the experimental curves are almost 
undist inguishable . 

Finally, more experiments at different fluid heights, surfactant concentrations, 
and temperatures were performed. Figure 4 gathers experimental values of 5/ 
inferred from these measurements as a function of bjth predicted by eq. (2) 
where Gq was extracted from linear rheology. This plot clearly shows that 
eq. (2) holds for all our experiments so that our interpretation in terms of 
standing elastic waves is relevant. 



3 Analytical and numerical predictions of the linear stability anal- 
ysis 



In this section, we compare our experimental results to the numerical approach 
first proposed by L. Tuckermann and K. Kumar (jsl) and adapted to viscoelastic 
fluids by S. Kumar (llSl ) by taking into account a complex, frequency dependent 
viscosity ri*\^\ 



3. 1 Numerical predictions for Oc and kc 



The numerical scheme relies on a linear stability analysis based on Floquet 
theory. Fourier expansions of the velocity field v and of the surface deformation 
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^ together with hnearization of the Navicr- Stokes equations lead to a, N x N 
eigenvalue problem, where corresponds to the cut-off in the Fourier expan- 
sions (we chose to use = 8, but iV = 4 already provides good estimates). 
The eigenvalues a(tc;o, k) correspond to marginal values of the acceleration i.e. 
for a given pulsation cuq — Stt/ of the excitation, surface modes with wave 
number k go unstable when a > a{u!o,k). For a fixed ujq, a = a{uJo,k) de- 
fines a set of so-called "resonance tongues" in the (a, k) plane. Each resonance 
tongue is associated to a different resonance frequency uj = nuJo/2 of the sur- 
face, where n is a positive integer. Odd (resp. even) values of n correspond to 
a subharmonic (resp. harmonic) response. 
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Fig. 5. (a) Linear rheology of a 3% wt. CPyCl-NaSal solution: experimental vis- 
coelastic moduli G' (•) and G" (o) versus shear oscillation frequency /. The dotted 
line and the continuous line respectively represent the best fits by a Maxwell fluid 
over / = 0-3 Hz (see table 1), and by an Oldroyd fluid over / = 0-40 Hz (see ta- 
ble 2). (b) Critical acceleration Qc and (c) critical wave number kc versus vibration 
frequency for a 3% wt. CPyCl-NaSal solution. The dotted (resp. continuous) lines 
are the result of the numerical calculation for the previous Maxwell (resp. Oldroyd) 
fluid. 

This linear stability analysis allows one to determine the critical acceleration Qc 
as the global minimum of the resonance tongues and to define the correspond- 
ing wave number as the critical wave number kc- In all the cases investigated 
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in the work, Oc and kc were found to lie on the first resonance tongue, so 
that the surface response is always subharmonic at u = ujo/2 = nf . The only 
severe constraint of the numerical method of ref. (jisl ) is that it requires the 
knowledge of f]*{uj) over a wide range of frequencies, typically O-Nf/2 ~ 2/ 
for an external forcing frequency /. Since / can reach 120 Hz in our experi- 
ments, rheological measurements of i]*{uj) are in principle required up to about 
250 Hz, which is far above the upper limit of our rheometer (about / = 40 Hz). 
Thus, above / ~ 20 Hz, extrapolations of ri*{uj) were used in the numerical 
calculation based on the two models discussed above: the Maxwell and the 
Oldroyd models. 



Figure 5(a) shows the best fit of the viscoelastic moduli of our 3% wt. CPyCl- 
NaSal solution by a Maxwell fiuid, which yields an accurate description of the 
experimental G' and G" up to / ~ 2 Hz. Feeding the numerical calculation 
of ref. (jlSl ) by this Maxwell model leads to the values of and kc shown as 
dotted lines in fig. 5(b) and (c). The first noticeable result is the presence of 
strong oscillations in Oc with the same periodicity as the experimental ones. 
Moreover the evolution of k^ is qualitatively similar to that of the experi- 
mental data. Still experiment and theory remain quantitatively different. The 
critical acceleration predicted numerically is smaller than the experimental Oc 
by about an order of magnitude and, for both Oc and kc, the amplitude of 
the oscillations predicted from the Maxwell model increases with frequency 
whereas it is clearly seen to decrease in the experiment. Since the vibration 
frequency is always greater than 20 Hz, these discrepancies are most probably 
linked to the poor fit of G" by the Maxwell model for / > 2 Hz. 



To account for viscoelasticity at higher frequencies, the rheological measure- 
ments were fit by an Oldroyd fiuid. As already pointed out, the Oldroyd model 
has two free parameters ti and T2 and leads to a better fit of the viscoelastic 
moduli for / > 2 Hz (see the continuous lines in fig. 5(a)). The linear stabil- 
ity analysis using the corresponding Oldroyd fiuid also yields a much better 
agreement between numerical and experimental Oc and kc. In particular the 
amplitude of the oscillations now decreases with increasing frequency as seen 
in the experiment. However, for / > 60 Hz, the critical acceleration found 
numerically remains smaller than the experimental one by about 20%. 

At this stage the lack of rheological data at frequencies larger than 40 Hz does 
not allow us to draw any definite conclusion on whether a better agreement 
between theory and experiment would be achieved by using a more elaborate 
rheological model at high frequencies. Experimental imperfections such as ab- 
sorption of shear waves at the bottom of the container or additional dissipation 
at the walls, which are not taken into account in the numerical calculation, 
may also account for the observed discrepancies. Another possibility is the fact 
that surface waves may affect the microstructure of the fluid in a nonlinear 
way (see sect. 5 below). Therefore we did not try to look further for more ap- 
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propriate rheological models but we rather focused on the original features of 
the Faraday instability in our strongly viscoelastic fluid, namely the presence 
of elastic waves across the sample. 



3.2 Calculation of the velocity field 



Interestingly the velocity field v of a vertically vibrated viscoelastic fluid at 
onset can be inferred from the numerical model of ref. (jlSl ). Indeed, if the 
velocity vector is expanded in terms of a Fourier series 



oo 

^ v„(a;,y,zje 

n=0 



i(n+o)a;ot 



(3) 



with a = (resp. a = 1/2) for a harmonic (resp. subharmonic) response, then 
ref. ( ,15i ) shows that the n"" Fourier mode of the vertical component w of the 
velocity reads 



Wn{x, y, z) = Vn{x, y, z) ■e^ = f{x, y)wn{z) (4) 
= /(x, y) (aie'^^" + a^e-^^-' + age^"^ + 046"'?"^) , (5) 



where 



+ \f{x,y) = ktf{x,y), (6) 



ri*[[n + ajujo) 
and kc is the critical wave number. 

The Oj's are deduced from the boundary conditions and depend on n, a, fee, t]* , 



u, and h (see ref. (1l5l ) for their full expressions). They are also proportional to 
the (n + a)* Fourier component of the surface wave amplitude However, 
since eq. (5) derives from a linear stability analysis, the velocity field is only 
known up to some multiplicative constant and a weakly nonlinear approach 
would be required to get a quantitative estimate of the velocity amplitude as 
a function of the distance from instability threshold. In the present work, we 
shall assume that our experiments are performed close enough to onset so as 
to allow us to use eq. (5) where the absolute value of the velocity will be left as 
a free parameter. Finally, the horizontal velocity ue^ + vey can be calculated 
from w and the incompressibility equation V ■ v = 0. 
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Fig. 6. (a) Numerical velocity field v(x, z, to = 0) for an Oldroyd fluid with 
rjo = 20 Pa.s, ri = 0.5 s, and T2 = 0.5 ms vibrated at / = 70 Hz. (b) Phase 
(j){x,z) and (c) amplitude u{x,z) of the x-component of the velocity v{x,z,t). 



Figure 6 presents the first subharmonic mode {n — 0,a — 1/2) oi the velocity 
field calculated from eq. (5) for an Oldroyd fluid. The velocity field consists of 
an array of counterrotating rolls whose number across the fluid height increases 
with increasing frequency. Of course, these counterrotating rolls are not steady 
rolls but oscillate in time at the frequency of the surface response. As shown by 
eq. (7), oscillations along the z direction arise from the term 036^"^ + 046"^"^ 
in eq. (5). Indeed, since qn is a complex number, it can be rewritten as qn — 
ikn + I /hi, where kn is the vertical wave number and In is the attenuation 
length for the {n + a)"' mode. 

Focusing on the case {n = 0,a = 1/2) and neglecting higher-order modes in 
eq. (5), it is straightforward to understand why couterrotating rolls should 
be detected in experiments on wormlike micelles but should remain unseen 
in Newtonian (or weakly viscoclastic) fluids. Indeed, with kc = 600 m^^ and 
ujq = 380 rad.s~^ in standard experimental conditions and for an Oldroyd fluid 
with typical parameters rjo = 20 Pa.s, ri = 0.5 s, and T2 = 0.5 ms, one gets 
Iq — 16 mm, which is greater than the vertical wavelength 27r/ko — 9 mm. 
On the other hand, for a Newtonian fluid of same absolute viscosity 77 = 
|?7*(a;o/2)| = 0.43 Pa.s, the calculation yields Iq = 1.6 mm ^ 27r/A;o = 36 mm. 
Note also that in the latter case Iq < h = 10 mm so that elastic waves are 
fully damped when they reach the bottom of the cell and no stationary elastic 
wave pattern can form across the fluid height. 



13 



Finally eq. (7) also allows one to recover eq. (2) in the case of an Oldroyd 
fluid and when only the first subharmonic mode is considered. In such a case, 
inserting eq. (1), n = 0, and a = 1/2 in eq. (7) leads to: 



*^ = *2 + >^ii^, (8) 

r^o 1 + ^t^T2 

where uj = ujq/2 is the pulsation of the surface response. Since in our exper- 
iments U!r2 ^ 1 <^ UTi (see table 2) and k'^rjo / pu^Ti <^ 1, the zero-order 
approximation of eq. (8) yields: 




As explained above the oscillations of Oc and kc are linked to the wave number 
go- More precisely one period of the oscillation of corresponds to an increase 
of ko by n/h (addition of a local maximum of v) i.e. an increase of u by 
-k^Gq/ ph^. For a subharmonic response this leads to 5f = \JGq/ pK^, i.e. 
eq. (2). Thus our simple interpretation in terms of standing elastic waves is 
consistent with the analytical expression (7) at least in a first approximation. 

Expanding go to first order, one gets : 

. ry ikl [G^ ry l + rirau;^ 
q^^iuj — - — . 10 

y Gq 2uj \I p y Go 2ri 



With go = iko + l/lo, the previous equation leads to 



k, = uJ^-^J^, (111 




lo = \— ,_, ' (12) 

Since kc increases roughly linearly with uj, eq. (2) becomes less accurate as 
frequency increases due to the corrective term in eq. (11). This explains why 
maxima of Oc do not correspond to integer values of f /6fth in fig. 3(a). More- 
over eq. (12) shows that the attenuation length Iq decreases with uj. At high 
frequencies oscillations should thus disappear. This accounts qualitatively for 
the decrease of the amplitude of the oscillations in ac and kc observed experi- 
mentally and numerically. 



So far we have shown that the linear stability analysis of ref. (Il5l ) predicts an 
oscillating behaviour for ac and kc- The calculation of the velocity field close 
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to instability threshold unveils the existence of counterrotating rolls. Such 
rolls, unseen in Newtonian fluids, are a clear signature of viscoelasticity in the 
Faraday experiment. Our next step is to provide direct experimental evidence 
for the existence of these "elastic rolls." 



4 Experimental study of the velocity field 



In this section we flrst describe the experimental apparatus used to access the 
velocity fleld in the Faraday experiment. Measurements are then presented 
and discussed in light of the above theoretical predictions. 



4.1 Particle imaging velocimetry setup 



In order to provide experimental evidence for the existence of large-scale coun- 
terrotating rolls in our viscoelastic fluid, particle imaging velocimetry (PIV) 
experiments were performed under vibrations. Since wormlike micelle solutions 
are transparent, the scattering properties of our system were enhanced by seed- 
ing the solutions with 0.1% wt. hollow glass spheres (Sphericel, Potters Indus- 
tries) of mean radius 11.7 fim and density 1.1. Thanks to the high zero-shear 
viscosity of the micellar solutions, sedimentation of the spheres is negligible 
over the ~8 hour duration of a Faraday experiment (a rough estimate yields a 
sedimentation velocity of about 0.03 /im.s"^). We also checked that the rheo- 
logical properties of the micellar system were not signiflcantly affected by the 
addition of 0.1% wt. scatterers i.e. remain very well described by an Oldroyd 
fluid. The best flt parameters obtained for the Oldroyd model on a 4% wt. 
CPyCl-NaSal solution seeded with 0.1% wt. scatterers are rjo = 24.3 Pa.s, 
Ti = 0.87 s, and T2 = 0.55 ms, very close to those of the corresponding solu- 
tion without scatterers (see table 2). 

In order to visualize the motion of the scatterers under vibrations, we use a 
parallelepipedic Plexiglas cell of length L = 180 mm (in the x direction), width 
/ = 4 mm (in the y direction), and depth h = 10 mm (in the z direction) under 
the brimful boundary condition. Thanks to the large aspect ratio L/l, a one- 
dimensional surface wave pattern is forced perpendicular to the x direction 



(1321) . The cell is lit in the (x, z) plane by a vertical laser sheet crossing the cell 
at I/O = 2 mm. Images of the {x, z) plane are recorded by a high-speed CCD 
camera (Mikroton MC1310) at a frequency ~ 10/ during 40 periods of the 
vertical oscillations. 

The flrst step in the image processing is to detect the vertical position of 
the cell by following a bright spot made by a laser at the bottom of the cell 
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5 10 15 20 25 

X (mm) 

Fig. 7. Picture of a 4% wt. CPyCl-NaSal solution seeded with hollow glass spheres 
recorded under vibrations at / = 75 Hz by the high-speed CCD camera. The arrows 
show the velocity field ue^ + wBz computed by PIV on two consecutive images. 
The white spots around x 2± correspond to the laser beam used for detecting the 
bottom of the cell. 

(around ~ 0,2; ^ —10 mm in fig. 7, see also fig. 11c). The displacement 
field of the scatterers and the two-dimensional velocity vector ue^ + w at 
a given point {x, z) are then computed in the moving reference frame of the 
vibrated cell by cross-correlating the intensity fields of two consecutive images 
over small square regions of size 0.25 mm^. Figure 7 shows a typical image 
of the micellar solution along with the velocity field obtained by PIV under 
vibrations. 

For the sake of clarity and since y is fixed to yo in all the PIV experiments, we 
shall drop the dependence on y of the various variables defined in the previous 
section whenever such dependence is not needed. Note also that the ±1 pixel 
uncertainty on the vertical position of the cell may remain too large compared 
to the vertical displacement of the fluid and in some cases leads to artifacts 
in the estimation of the vertical component w of the velocity. Therefore, in 
the following, we focus on the velocity component u along the x direction. In 
order to compare experimental results to numerical calculations, we extract 
the mode oscillating at ujo/2 from u{x,z,t) i.e. Uq{x,z) as defined by eq. (3) 
with a — 1/2. 



4-2 Experimental velocity measurements compared to numerical predictions 



Shown in fig. 8 are the amplitude u{x,z) and phase (f){x,z) of Uo{x,z). It is 
clearly seen that the velocity does not decrease exponentially with z as would 
be expected for a Newtonian fluid, but rather presents local maxima consistent 
with the presence of counterrotating rolls (see fig. 6(b) and (c) for a qualitative 
comparison with numerical results). 

In the case of our one-dimensional pattern with large aspect ratio, one may re- 
cover the horizontal velocity u from the numerical calculation of w by assuming 
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Fig. 8. (a) Phase (j){x,z) and (b) amplitude u{x,z) of uo{x,z) measured by PIV in 
a 4% wt. CPyCl-NaSal solution vibrated at / = 75 Hz. (c) < z) > and (d) 
< u{x, z) > averaged over x = 6-9 mm. The continuous lines are the predictions 
for (j) and u given by equations (5) and (13) with n = and a = 1/2 and using 
an Oldroyd fluid with tjq = 24.3 Pa.s, ri = 0.87 s, and T2 = 0.55 ms. Both the 
amplitude of uo and the reference phase are free parameters. 



L ^ I. Indeed under this assumption and the no slip condition at the cell bor- 
der (v = for 1/ = 0, /) one gets v = and w{x, y, z) = cos{kxx) cos{kyy)w{z), 
with ky = tt/I and k^ = k^ — ky. This is equivalent to having a half mode 
along the y direction. Next, incompressibility leads to: 



Un{x, y, z) = -Y sm{kxx) cos{kyy) ^^(^) , (13) 



where Wn{z) is given by eq. (5). 

Figures 8(c) and 8(d) compare the experimental results (averaged over x = 6- 
9 mm) to the numerical prediction of eq. (13) for [n = 0,a = 1/2) corre- 
sponding to the best fit of the hnear rheological measurements by an Oldroyd 
fluid. The number of rolls as well as their positions in the cell are very well 
accounted for by the numerical calculation. Wc believe the discrepancies near 
the surface, especially the large experimental values of Uq, to be linked to the 
non-zero amplitude of the surface waves that prevents an accurate estima- 
tion of the velocity field. 
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5 Birefringence experiments 



Measurements of Oc and fee, experimental characterization of tlie velocity field 
close to onset, and comparison with numerical predictions have confirmed 
the interplay of parametrically excited surface waves and elastic waves in a 
strongly viscoelastic fluid under vertical vibrations. Thus the presence of a 
microstructure, namely an entangled network of wormlike micelles, deeply af- 
fects the Faraday instability. We now address the question raised in ref. ( jisl ) 
of whether the instability modifles the microstructure. In particular, in the 
nonlinear regime, semi-dilute solutions of wormlike micelles are known to dis- 
entangle and align under shear leading to a shear-induced isotropic-nematic 



transition (Il9l : l33l : |3J). Therefore one may wonder whether the flow induced 
by Faraday waves is strong enough to induce such a transition, which could 
provide an alternative explanation for the discrepancy between experiments 
and linear stability analysis observed in sect. 3.1. 

Flow-induced alignment have already been much studied in semi-dilute worm- 
like micelle solutions using birefringence S; S; M, S H 0; iH) . Following 



these previous works and in order to investigate the deformation of the mi- 
crostructure under vertical vibrations, we performed birefringence experiments 
close to onset and further above the instability threshold. The experimental 
setup is modifled as follows. A parallelepipedic glass cell of length L = 72 mm 
and width / = 6 mm is flUed with the mi cellar solution to a height h = 7.1 mm. 
This cell is set between two polarizers and lit from backward by an extended 
white light source. The high-speed CCD camera used previously for PIV mea- 
surements captures the transmitted light at about 1000 fps. For given vibration 
frequency / and acceleration a, two measurements can be performed depend- 
ing on the conflguration of the polarizers. When the polarizers are aligned, 
the contour of the surface wave is easily accessed by image processing. This 
yields the surface deformation ^{x, t) from which the amplitude of the surface 
wave ^0 may be extracted. When the polarizers are crossed, the birefringence 
intensity I{x,z,t) is recorded and analyzed. 



5.1 First observations and qualitative discussion 



Figure 9 shows a typical birefringence pattern obtained at a given time to 
between crossed polarizers on an 8% wt. CPyCl-NaSal solution slightly above 
the onset of Faraday waves. Evenly spaced bright spots fill the whole cell with 
the same horizontal periodicity as the surface wave. Such bright spots are 
observed below antinodes of the stationary surface wave. The birefringence 
intensity I{x,z,t) oscillates in phase with ^{x,t) so that bright spots alterna- 
tively turn to black when the surface deformation goes through zero. Regions 
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Fig. 9. Typical birefringence pattern I{x,z,to) of an 8% wt. CPyCl-NaSal solution 
vibrated at / = 80 Hz. 




45 50 55 60 65 70 75 80 
f(Hz) 



Fig. 10. (a) Number of birefringent spots Ugpot along the vertical direction and (b) 
wave number kx versus vibration frequency for a 5% wt. CPyCl-NaSal solution. 

in the bulk located below nodes of the surface wave remain dark at all times. 

Based on the previous analysis of the velocity field, the interpretation of 
the birefringence pattern is rather straightforward. Indeed, from the veloc- 
ity fields of figs. 6 or 7, the two-dimensional rate of deformation tensor 7 
can easily be estimated using 7,jj = diVj, with i and j being x, y, or z and 
{vx,Vy,Vz) = {u,v,w). In the {x,z) plane, 7 (not shown) is composed of local- 
ized compression and stretching zones that oscillate in time. Compression and 
stretching result in an oscillatory deformation of the micellar network, which 
leads to the observed birefringent spots. 

Moreover, as shown in fig. 10(a), the number Ugpot of bright spots along the 
vertical direction increases with increasing frequency. Figure 10(b) proves that 
drops of kx are linked to the discretization of the number of spots along the 
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z direction. This is a direct consequence of the hnk between the vertical wave 
number and (see eq. (7)): each time the number of spots increases by one 
(or equivalently each time a new set of counterrotating rolls fits in the cell 
height), the vertical wave number increases (by 7r/2/i) and the surface wave 
number abruptly decreases. 

Finally, if one defines an equivalent shear rate near the surface as 7 = 



(Il5l ). a rough estimate with ,^0 = 30 /im (the minimum amplitude that our 
setup can detect) yields 7 ~ 3 s~^ Since such a shear rate is of the same order 
of magnitude as the critical shear rate 7c = 2.6/ri for micelle disentanglement 
and alignment (34^), nonlinear rheology may come into play and we should 
check for the possibility of an alignment transition of the micelles under ver- 
tical vibrations. The next sections are devoted to a detailed analysis of the 
birefringence intensity and to more measurements above onset. 



5.2 Birefringence measurements above onset 



For an optically anisotropic medium composed of elongated and oriented mi- 
celles, the transmitted birefringence intensity takes the following general ex- 



pression (1371 ): 



/(x, z, t, 6) = I' sin^ ( ^i^lM ) sin^ 2{x{x, z, t) - 9) , (14) 



where the retardation S reads 
2ttI 

6{x, z,t) = ——An{x, z,t), (15) 
A 



/' is the incident intensity, A ~ 500 nm is the light wavelength, and An{x, z, t) 
is the birefringence intensity, xi^y ^5 1) is defined as the angle between some 
reference direction and the average micelle direction. 6 is the angle between 
the polarization of the incident beam and the reference direction. 

To address the issue of whether or not the nonlinear deformation regime is 
reached and the micelles disentangle, we shall make use of the so-called "stress- 
optical rule" according to which the quantity An sin 2x is proportional to the 



shear stress a (1371 : l38l ). This rule was shown to hold at low shear stresses in 
CTAC wormlike micelle solutions for various salt (NaSal) concentrations (jssh . 
Although to our knowledge no experimental validation of the stress-optical 
rule is available for our CPyCl-NaSal in brine system, we may assume that it 
is general enough to also hold in our experiments. 
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Fig. 11. (a) Surface wave amplitude ^(xo,t), (b) birefringence intensity I{xq, ZQ,t), 
and (c) position of the bottom of the ceh b{t) for an 5% wt. CPyCl-NaSal solution 
vibrated at / = 65 Hz. In (b) xq = 7.5 mm and zq = —2.1 mm were chosen to 
coincide with a birefringence spot. 

As long as the network structure remains intact and the micelles remain 
entangled, the deformation is elastic, i.e. a oc 7, so that one should have 
Ansin2x oc 7, where 7 stands for the deformation. However the stress-optical 
rule no longer applies when the micelles disentangle and the network is broken 



(137|;138D. 



In the case 5 and for ^ = 0, which is always verified for wormlike micelles 



(1371) ■ eq. (14) and the stress-optical rule lead to 

/(x, 2;, t) oc An(x, tY sin^ 2x(x, t) oc 7(x, t)^ . (16) 



Focusing on the first subharmonic mode, eq. (16) shows that I{x,z,t) should 
oscillate at the forcing pulsation uq with an amplitude Io{x,z) oc 70(0;, z)^, 
where 7o(a;, z) denotes the amplitude of the first subharmonic mode of 7(x, z, t). 
Experimentally 70 is varied by increasing the driving acceleration a above the 
critical acceleration Oc. Indeed the deformation at the surface ■jo{x, z = ^oi^)) 
is directly proportional to the surface wave amplitude through 70 = fcc^o- Thus 
increasing a leads to larger surface wave amplitudes and to larger deforma- 
tions, which allows one to test eq. (16). 
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Fig. 12. (a) Surface wave amplitude versus acceleration a for an 8% wt. CPyCl — 
NaSal solution vibrated at / = 80 Hz. (b) versus a. The dotted line is the best 
fit by eq. (19) over a = 88.3-90.1 m.s~^ and the continuous line is the best fit by 
eq. (20) over a = 88.3-97 m.s-^. 



More precisely, we first extract the position of the surface C,{x,t) from images 
recorded with aUgned polarizers. The subharmonic amplitude of the deforma- 
tion C,o{x) is then computed and fitted to ^0(2;) = ^o| cos^k^x)] (see fig 11(a)), 
which yields the amplitude of the deformation 70 = fcc^o at the surface. Since 
the dependence of 70(^,2;) on z is decoupled from the dependence on x (see 
also eq. (4)), one has 7o(x, z) = 7o| cos{kxx)\f{z), where f{z) does not need 
to be specified in order to test the scaling law (16). 



Figure 12 shows the evolution of the surface wave amplitude ^0 with the driving 
acceleration a above onset. In the case of the Faraday instability, it is well 
known that ^0 should obey the following amplitude equation (1321 ) : 



dt 



■So 



(17) 



where Tg is the characteristic growth time of the instability, e = a/ac — 1 is 
the reduced acceleration, and I3 and I5 are characteristic lengths linked to the 
dissipation process. At steady state eq. (17) becomes: 



f3 t5 

^0 so 
~ 72" ~ 74 + ••• 



0. 



22 



If the expansion is truncated to third order, one finds the usual expression: 



^0 = ^av^- 



(19) 



As seen in fig. 12, eq. (19) accounts for the experimental data over a rather 
restricted range of accelerations and the term should be included (jsl). In 
this case, solving eq. (18) leads to: 



which yields a much better description of the .^o measurements. 

Once ^0 (hence 70) is known, the birefringence intensity I{x,z,t) is measured 
from images obtained between crossed polarizers at ^ = so that the intensity 
at the bright spots is maximized. Figure 11(b) shows I{xo,ZQ,t) recorded at 
the position of a bright spot. Contrary to the prediction of eq. (16), I{xo, zq, t) 
is not perfectly harmonic but contains a significant subharmonic component. 
This is most probably due to the finite amplitude of the surface waves which 
induce a small displacement of the position of the birefringence spots at fre- 
quency f/2. To decrease the noise in the determination of /q, I{xo,zo,t) is 
averaged over two horizontally adjacent bright spots and Jq is defined as the 
harmonic component of this mean intensity. 

Figure 13 presents the evolution of Iq as a function of 70 for two different 
concentrations and driving frequencies. We find that Iq oc 7q up to a = l.loc. 
Beyond this point the surface becomes chaotic. Our conclusion is that the 
stress-optical rule and eq. (16) hold and that the micellar network does not 
break, at least for a < l.la^. In particular micelles remain entangled close 
to the instability threshold and nonlinear rheology cannot be invoked to ex- 
plain deviations from linear stability analysis. Thus, even if the amplitude 
of the equivalent shear rate falls into the nonlinear regime, it appears that 
the micellar network remains intact under vibrations. We suggest that this 
is linked to the short period of the forcing (1// ~ 0.02 s) compared to the 
characteristic time of our fiuid (ri ~ 0.5-1 s). Rather than the amplitude 
of the equivalent shear rate, it would probably be more relevant to compare 
the deformation ^o^c to the critical deformation ri7c for disentanglement and 
alignment. A good approximation is ri7c — 2.6 while ^o^c — 0.2 for a = l.loc. 
Since ^o^c -C Tijc, the fiuid deformation over one period of the forcing is too 
small to break the network and induce the alignment transition. 




(20) 
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Fig. 13. Square root of the transmitted intensity Iq versus the deformation 70 for 
a 5% wt. CPyCl-NaSal solution vibrated at / = 65 Hz (o) and for a 8% wt. 
CPyCl-NaSal solution vibrated at / = 80 Hz (•). The continuous hne is the first 
diagonal, (a) Linear scales, (b) Logarithmic scales. 

Conclusion 



In this paper we have presented new results for the Faraday instabihty in 
wormlike micelles. We have shown that the strong viscoelasticity of these 
solutions gives rise to original features characterized by oscillations of the 
critical acceleration and wave number. We clearly linked these oscillations to 
the presence of standing elastic waves across the cell. Velocimetry experiments 
have provided evidence that such elastic waves generate count errot at ing rolls 
so that the velocity field does not decrease exponentially with depth. Linear 
stability analysis was shown to predict and fit both these phenomena. Finally 
birefringence measurements have shown that the deformation generated at the 
interface tends to orient the micellar network but is not large enough to induce 
disentanglement . 

We believe that the present experiments may contain worthwile implications 
for rheological studies of strongly viscoelastic fluids. As noted above and as 
already mentioned in ref. (fiol ). the Faraday instability constitutes an easy way 
to probe a complex fluid in an oscillating elongational flow at frequencies above 
30 Hz. One may thus wonder whether rheological information can be gained 
from measurements under vertical vibrations and how the Faraday experiment 
may be used as a "high" frequency rheometer. For instance, we tried to infer 
the rheological properties of our micellar solutions by fitting our velocity field 
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measurements to equations (3-7) without assuming any rheological model 



(l3ll). Although the results did not prove conclusive enough (probably due to 
boundary effects as mentioned at the end of sect. 3.1), this may still constitute 
an interesting direction for future research. Finally another perspective for this 
work would be to study complex fluids with shorter characteristic times, e.g. 
liquid crystals, in which the surface waves may induce microstructural changes. 

The authors wish to thank P. Lettinga and B. Pouligny for fruitful discus- 
sions. S. Lerouge is thanked for advice on the birefringence experiments and 
P. Snabre for technical help on the PIV measurements. 



References 

[I] M. Faraday, On the forms and states assumed by fluids in contact with 
vibrating elastic surfaces, Philos. Trans. R. Soc. Lond. 52 (1831) 319. 

[2] Lord Rayleigh, On the crispations of fluid resting upon a vibrating sup- 
port, Phil. Mag. 16 (1884) 50. 

[3] T. B. Benjamin and F. Ursell, The stability of the plane free surface of a 
liquid in a vertical periodic motion, Proc. R. Soc. London A 225 (1954) 
505. 

[4] P. Chen and H. Vinals, Pattern selection in Faraday waves, Phys. Rev. 
Lett. 79 (1997) 2670. 

[5] W. S. Edwards and S. Fauve, Patterns and quasi-patterns in the Fara- 
day experiment, J. Fluid Mech. 278 (1994) 123; J. Bechhoefer, V. Ego, 
S. Manneville, and B. Johnson, An experimental study of the onset of 
parametrically pumped surface waves in viscous fluids, J. Fluid Mech. 
288 (1995) 325; E. A. Cerda and C. T. Tirapegui, Faraday's instability 
in viscous fluid, J. Flmd Mech. 368 (1998) 195. 

[6] A. Kudrolli and J. P. Gollub, Superlattice patterns in surface waves, 
Physica D 97 (1996) 133. 

[7] P. Manneville, Dissipative Structures and Weak Turbulence (Academic 
Press, San Diego, 1990) 

[8] K. Kumar and L. S. Tuckerman, Parametric instability of the interface 
between two fluids, J. Fluid Mech. 279 (1994) 49. 

[9] A. Wernet, C. Wagner, D. Papathanassiou, H. W. Miiller, and K. Knorr, 
Amplitude measurements of Faraday waves, Phys. Rev. E 63 (2001) 
036305. 

[10] F. Raynal, S. Kumar, and S. Fauve, Faraday instability with a polymer 
solution, Eur. Phys. J. B 9 (1999) 175. 

[II] C. Wagner, H. W. Miiller, and K. Knorr, Faraday waves on a viscoelastic 
liquid, Phys. Rev. Lett. 83 (1999) 308. 

[12] F. S. Merkt, R. D. Deegan, D. I. Goldman, E. C. Rericha, and H. L. Swin- 
ney. Persistent holes in a fluid, Phys. Rev. Lett. 92 (2004) 184501. 



25 



p. Huber, V. P. Soprunyuk, J. P. Embs, C. Wagner, M. Deutsch, and 
S. Kumar, Faraday instability in a surface- frozen liquid, Phys. Rev. Lett. 
94 (2005) 184504. 

A. V. Kityk and C Wagner, Delay of disorder by diluted polymers, Eu- 
rophys. Lett. 75 (2006) 441. 

S. Kumar, Parametric driven surface waves in viscoelastic liquids, Phys. 
Fluids 11 (1999) 1970; Vibration-induced interfacial instabilities in vis- 
coelastic fluids, Phys. Rev. E 65 (2002) 026305. 

H. W. Miiller and W. Zimmermann, Faraday instability in a linear vis- 
coelastic fluid, Europhys. Lett. 45 (1999) 169. 

O. Lioubashevski, Y. Hamiel, A. Agnon, Z. Reches, and J. Fineberg, 
Oscillons and propagating solitary waves in a vertically vibrated colloidal 
suspension, Phys. Rev. Lett. 83 (1999) 3190. 

P. Ballesta and S. Manneville, Signature of elasticity in the Faraday in- 
stability, Phys. Rev. E 71 (2005) 026308. 

R. G. Larson, The Structure and Rheology of Complex Fluids (Oxford 
University Press, Oxford, 1999). 

C.-H. Liu and D. J. Pine, Shear-induced gelation and fracture in micellar 
solutions Phys. Rev. Lett. 77 (1996) 2121. 

J. F. Berret, G. Porte, and J. P. Decruppe, Inhomogeneous shear flows 
of wormlike micelles: a master dynamic phase diagram, Phys. Rev. E 55 
(1997) 1668. 

J.-F. Berret, D. C. Roux, and G. Porte, Isotropic-to-nematic transition 
in wormlike micelles under shear, J. Phys. II France 4 (1994) 1261. 
J.-F. Berret, J. Appell, and G. Porte, Linear rheology of entangled worm- 
like micelles, Langmuir 9 (1993) 2851. 

M. E. Gates, Reptation of living polymers: dynamics of entangled 
polymers in the presence of reversible chain-scission reactions. Macro- 
molecules 20 (1987) 2289. 

H. Rehage and H. Hoffmann, Rheological properties of viscoelastic sur- 
factant systems, J. Phys. Chem. 92 (1988) 4712. 

R. B. Bird, R. C. Armstrong, and O. Hassager, Dynamics of Polymeric 
Liquids (Cambridge University Press, Cambridge, 1987). 
O. Manero, F. Bautista, J. F. A. Soltero, and J. E. Puig, Dynamics of 
worm-like micelles: the Gox-Merz rule, J. Non-Newtonian Fluid Mech. 
106 (2002) 1. 

B. Yesilata, C. Clasen, and G. H. McKinley, Nonlinear shear and exten- 
sional flow dynamics of wormlike surfactant solutions J. Non-Newtonian 
Fluid Mech. 133 (2006) 73. 

I. A. Kadoma, G. Ylitalo, J. W. van Egmond, Structural transitions in 
wormhke micelles, Rheol. Acta 36 (1997) 1. 

P. Fischer and H. Rehage, Rheological master curves of viscoelastic surfac- 
tant solutions by varying the solvent viscosity and temperature, Langmuir 
13 (1997) 7012. 

P. Ballesta, Instabilite de Faraday dans les fluides complexes (Ph. D. 



26 



Thesis, Univcrsite Bordeaux 1, unpublished, 2006). 
[32] S. Douady, Experimental study of the Faraday instability, J. Fluid. Mech. 
221 (1990) 383. 

[33] H. Rehage and H. Hoffmann, Viscoelastic surfactant solutions: Model 
systems for rheological research, Mol. Phys. 74 (1991) 933. 

[34] N. A. Spenley, M. E. Gates, and T. C. B. McLeish, Nonlinear rheology 
of wormlike micelles, Phys. Rev. Lett. 71 (1993) 939. 

[35] Y. Hu, S. Q. Wang, and A. M. Jamieson, Rheological and flow bire- 
fringence studies of a shear-thickening complex fluid-A surfactant model 
system, J. Rheol 37 (1993) 531. 

[36] E. Cappelaere, J. F. Berret, J. P. Decruppe, R. Cressely, and P. Lindner, 
Rheology, birefringence, and small-angle neutron scattering in a charged 
micellar system: evidence of a shear-induced phase transition, Phys. Rev. 
E 56 (1997) 1869. 

[37] S. Lerouge, J. -P. Decruppe, J.-F. Berret, Correlations between rheological 

and optical properties of a micellar solution under shear banding flow, 

Langmmr 16 (2000) 6464. 
[38] J. P. Decruppe and A. Ponton, Flow birefringence, stress optical rule and 

rheology of four micellar solutions with the same low shear viscosity, Eur. 

Phys. J. E 10 (2003) 201. 
[39] S. Lerouge, J.-P. Decruppe, and P. D. Olmsted, Birefringence banding in 

a micellar solution or the complexity of heterogeneous flows, Langmuir 

20 (2004) 11355. 

[40] S. Chen and J. P. Rothstein, Flow of a wormlike micelle solution past a 
falling sphere, J. Non-Newtonian Fluid Mech. 116 (2004) 205. 

[41] B. A. Schubert, N. J. Wagner, and E. W. Kaler, Shear-induced phase 
separation in solutions of wormlike micelles, Langmuir 20 (2004) 3564. 



27 



